Powerful Carnot cycle and attainability of the nicixinial efficiency 
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We want to understand whether and to which extent the maximal (Carnot) efficiency for heat 
engines can be reached at a finite power. To this end we generalize the Carnot cycle so that it is not 
restricted to slow processes. We show that for realistic (i.e. not purposefully-designed) engine-bath 
interactions, the optimal engine performing the generalized Carnot cycle has a long cycle time and 
hence vanishing power. This aspect is shown to relate to the theory of computational complexity. A 
physical manifestation of the same effect is the Levinthal's paradox in the protein folding problem. 
The resolution of this paradox for realistic proteins allows to construct engines that can extract at 
a finite power 50 % of the maximally possible work reaching 90 % of the maximal efficiency. For 
purposefully-designed engine-bath interactions, the Carnot efficiency is achievable at a large power. 



Reciprocating heat engines extract work operating 
cyclically between two thermal baths at temperatures Ti 
and T2 (Ti > T2) [l|. They have two basic characteris- 
tics: (i) efficiency, 77 = W/Qi, is the work W extracted 
per cycle divided by the heat input Qi from the high- 
temperature bath, (ii) Power W/t, where r is the cycle 
duration. Both these quantities have to be large for a 
good engine: if 77 is small, lot of energy is wasted; if the 
power is smalL no sizable work is delivered over a rea- 
sonable time [1|. 

The second law establishes the Carnot efficiency 77c = 
1 — 2k as an upper bound for 77 [l|. The Carnot cy- 
cle reaches the bounding value 77c in the (useless) limit, 
where the power goes to zero [l| . Conversely, realistic en- 
gines are not efficient, since they have to be powerful, e.g. 
the efficiency of Diesel engines amounts to 35-40 % of the 
maximal value. T^his power- efficiency dilemma TuotwaXcd 
a search for the efficiency that would generally charac- 
terize the maximal power regime. One candidate for 
this is the Curzon-Ahlborn efficiency 77CA = 1 ^ yT^jTi 
[2], which is however crucially tied to the linear regime 
Ti « T2 [si |j] . Beyond this regime rycA is a lower bound 
of 77 for a class of model engines [Sl . Several recent mod- 
els for the efficiency at the maximal power overcome 7ycA 
with ,,* = ^ @. 

As argued in [5|,lZ|,|8|, the maximal power regime allows 
for the Carnot efficiency, at least for certain models. But 
it is currently an open question whether the maximal ef- 
ficiency is attained under realistic conditions (see e.g. [9| 
in connection with [7| ) , and how to characterize the very 
"realism" of those conditions. Even more generally: what 
is the origin of the power-efficiency dilemma? We answer 
them by analyzing a generalized Carnot cycle, which in 
contrast to the original Carnot cycle is not restricted to 
slow processes. We now summarize our answers. 

(1) When the iV-particle engine operates at the max- 
imal work extracted per cycle, its efficiency reaches the 
Carnot bound 77c for iV ;:^ 1, while the cycle time is 
given by the relaxation time of the engine. The maximal 



work and the Carnot efficiency are achieved due to the 
flat energy landscape of the engine. For realistic engine- 
bath interactions this energy landscape leads to very long 
[©(e^)] relaxation time nullifying the power. By realistic 
we mean interactions that are independent from the en- 
gine Hamiltonian. If we assume a proper tuning between 
engine-bath interaction and the engine Hamiltonian, the 
maximal efficiency is achievable at a large power 0{yN). 

(2) The relaxation of the optimal engine under realistic 
interactions relates to an important problem of search- 
ing an unstructured database for a marked item, where 
each energy level refers to a database item. This task 
is computationally complex, i.e. even the most powerful 
quantum algorithms resolve it in 0{e^''^) time-steps |11| . 
Hence the power- efficiency dilemma relates to computa- 
tional complexity. The same effect can be reformulated 
as the Levinthal's paradox of the protein folding prob- 
lem: if the majority of unfolded states of a protein are 
assumed to have the same (free) energy, the folding time 
is very long J12l |. 

(3) A scenario of resolving the Levinthal's paradox pro- 
posed in protein science shows the way of constructing 
sub-optimal engines that operate at a reasonably large 
values of work, power and efficiency. These sub-optimal 
engines function as model proteins. Sacrificing some 50- 
60% of the maximal possible work leads to a reasonable 
relaxation times with the efficiency that achieves some 
90 % of the maximal (Carnot) efficiency. 

Carnot cycle and its generalization. Recall that 
the Carnot cycle consists of four slow, consecutive pieces 
[l|: thermally-isolated - isothermal - thermally-isolated 
- isothermal. Four times slow brings in the vanish- 
ing power stressed above. Since the overall process is 
a quasi-equilibrium one, the external fields that extract 
work from the engine act on it during all the four pieces. 
One deduces for the isothermal parts: Qi = TiAS and 
Q2 = T2AS, where Qi (Q2) is the heat taken from (put 
into) the Ti-bath (T2-bath), and AS* > is the entropy 
change. Since the work extracted is W ^ Qi — Q2, rj 



equals to its maximal value 77c = 1 — p- [l|- 

We keep the two isothermal and two thermally isolated 
pieces of the original Carnot cycle, but do not force them 
to be slow. In addition, the external fields will act only 
during the thermally isolated stages. Isothermal pieces 
amount to free relaxation; see Fig. 1. Due to this point, 
we shall be able to analyze the engine functioning from 
the energy conservation. We analyze the engine via quan- 
tum mechanics on a finite Hilbert space, because the in- 
herently discrete character of quantum mechanics makes 
possible the optimization of the cycle (i.e. it reduces the 
problem to a combinatorial optimization). The final re- 
sults arc interpreted classically and can be also obtained 
by discretizing the Hamiltonian classical dynamics over 
phase-space cells. 

0. The engine E with the Hamiltonian Hi starts in 
an equilibrium state at temperature Ti described by the 
density matrix 



p(0) = pi = e-'^i^7(tre-'^i^0, /?! = V^] 



-piHi 



(1) 



1. Between times and r, E undergoes a thermally iso- 
lated process with a time-dependent Hamiltonian Hi2{t) 
and the unitary evolution p(t) = Ui2p{Q)U{2'- 

Hr2{0) = Hi, Hi2{r)=H2, Uu ^Te'^ ^0 '''"''^'\{2) 

where T means chronological ordering. The work taken 
out of E is determined by energy conservation [see [l| and 
section I of the Supplementary Material] 



Wi = tiiHipi - H2Ur2PiUl2 



(3) 



2. Then E is attached to the T2-bath and after relax- 
ation time Tf its density matrix becomes [see Fig. 1] 



p{r + r,) = p2 = e-^^^V(tre-^=^=). 



(4) 



The heat that came to E from the r2-bath is 

O2 = tr[i72P2 - H2Ui2PlUl2]- (5) 

3. E undergoes another thermally isolated process 

iJ2i(0) - H2, H2i{t) ^ Hi, U21 - re-'^o<i-HM^)^{6) 

completing the cycle with respect to the Hamiltonian. 
The work taken out of E reads 



W2=tY[H2P2-HiU2lP2Ull]- 



(7) 



4. Finally, E is attached to the Ti-bath (Ti > T2) and 
relaxes to pi thereby completing the cycle; see ([T]). The 
heat that came to E from the Ti-bath is 



0i=tr[i7ipi-i7iC/2iP2t/Ji]- 



(8) 



To stress the differences with the original Carnot cycle: 
(i) the cycle time 2{t + Ti) need not be much larger than 



the relaxation time Tf. (ii) The cycle is out of equilib- 
rium. (Hi) The work source and the bath never act simul- 
taneously; either one acts or another; see Fig. 1. Hence 
heat and work are deduced from the energy conservation. 
The essential assumption of the cycle is that the 
system-bath coupling is weak; thus we did not count the 
work necessary for switching the coupling on and off. 
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FIG. 1: Generalized Carnot cycle. Full (dashed) arrows refer 
to isothermal (thermally isolated) pieces. Thermal baths are 
denoted by T\ and T2. Four boxes carry notations for the 
state and Hamiltonian; see ([THH]). The equilibrium states are 
pi ex e~"^''^^ and p2 oc e~"^''^^. Also, pi2 = U^piUl^, 
P21 = U2ip2Ul-^ are the states before isothermal pieces. 

Maximization of work. We maximize the full ex- 
tracted work W — Wi + W2 over Hi, H2, C/12, U21 for 
fixed Ti > T2 and a fixed number n -f 1 of energy levels of 
E. The lowest energies of Hi and H2 can be set to zero. 
Then W = Wi+W2 reads from (H ^ 
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where the eigenvalues and eigenvectors of Ha are 



,n+l 



and where, given a,"f — 1,2 and fc, I — 1, ..., n + 1, 
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Cli arc doubly stochastic matrices: Cj.^ > 0, 

Eltlcf^ = ELY4r' = 1- Every doubly stochas- 
tic matrix Cki can be represented as a convex sum of 
permutation matrices (Birkhoff's theorem) [1J|: Cki = 

T,s>^s'n^kl where A^ > 0, J^s^s = 1, and where hI-'I 
permutes the components of any vector on which it acts. 
Using Birkhoff's theorem we maximize ly in ([51 fTU)) . 
The optimal Cj.j and C^^ amount to permutation ma- 
trices, since they enter linearly into W. Hence they are 
unit matrices; see (fTUl) . In contrast to the original Carnot 
cycle, the optimal thermally isolated processes should 



proceed very quickly, so that the Hamiltonian changes, 
but the density matrix does not; sec (l3l[7]). The work 



>n+l 



W = W, + W2 = Y.u^y^ -P')^'i^ 4'0, (14) 

is to be stih maximized over {Cfc }fc^2 ^^^ {^i }fci2; 
see ([T^. W is symmetric with respect to permutations 

within {el^}ltl and within {e^^^}ltl- We checked nu- 
mericaUy that this symmetry is not broken and hence 
the maximum of W is reached for 



ejf' , a = 1,2, 



(15) 



i.e. ah excited levels have the same energy. With new 
variables e"^"*^ = Uq, we write the maximal work as 

(Ti In J- - T2 In Miui - U2)n 

max[Ml,W2] = r— tt—- , (16) 

[1 + nui\[l + nu2\ 
where ui and U2 arc found from 



9„iH/'max[wi,W2] = ^^^ Wniax ["l , "2] == 0. 



(17) 



Ml and U2 depend on T2/T1 and on n. Noting (|5]) and 
the result before ([T4)) we obtain Qi = tr{Hi{pi — P2)) 
for the heat obtained from the high-temperature bath. 
Using (jlSp we get from rj = W/Qi and from P^ : 



7? = 1 - [^2 lnu2] /[Ti In Ml] . 



(18) 



Note from dH]) that W,„ax[wi,'«2] > and T2 < Ti imply 
1 > £[2] /ell] > r2/ri. Hence ^ implies rj < rjc ^ 
1 — T2/T1, as expected. 

Thus the thermally isolated pieces of the cycle consist, 
respectively, of sudden changes e'^l -^ el^l and e'^' — >■ e'^l. 
This is an advantage, since sudden changes are automat- 
ically thermally isolated (no time to interact with the 
baths). In contrast, a slow thermally isolated change is 
difficult to realize, which is one of the reasons for the 
original Carnot cycle being practically useless [l| . 

Both Wniax[wii "2] and 77 increase with n. For ln[n] 3> 1 
we get asymptotically from (|17p : 



(l-6l)lnH 

Ml = , U,2 = 



where 6 = T2/T1. This produces 



7iln[n](l-6')' 



(19) 



W„,ax[iii, W2] = iT2 - Ti) ln[n] - 0(l/ln[n]), (20) 
,j^ijc-0{l/Hn]), Tjc^l-T2/Ti. (21) 

The maximal work Wmax[Mi, ^2] scales as ln[n], since this 
is the "effective number of particles" for the engine. In 
the macroscopic limit ln[n] 3> 1, the efficiency converges 
to its maximal value ijc = 1 — T2/T1; see (piT) . Recall 
from P^ that nua is the overall probability of the excited 
levels. Hence for ln[n] ^ 1 the ground state (excited 
levels) dominates in p2 (pi); see ([TI 



The cycle time amounts to two times the relaxation 
time Tf of the system with spectrum (jlSp and energy gap 
e ^ ln[7i]; see (|15[[TO)) . (Recall that the thermally isolated 
stages of the cycle are very quick.) The magnitude of t,- 
essentially depends on the scenario of relaxation. 

First (specific) scenario. We can assume that the 
Hamiltonian (jTTl [T5|) of the heat engine is known. Then 
there exist system-bath interaction scenarios that lead 
to Tr = 0{\/\n[n]); section IV of the Supplementary Ma- 
terial. Hence for this type of relaxation the Carnot ef- 
ficiency is achievable at a large power 0{^/]n\n\) ^ 1; 
see (pn|) . However, in these scenarios the system-bath 
interaction Hamiltonian (that governs the relaxation) is 
special: it depends on the engine Hamiltonian ([TTl [T5|) . 

Second (realistic) scenario. Assuming that the system- 
bath interaction does not depend on the Hamiltonian (|1H 
lisp , we can estimate Tj- within the weak-coupling, Markov 
master equation approach that leads to r^ = 0{n); sec- 
tion II of the Supplementary Material. For a qualitative 
understanding of this situation, consider the relaxation 
as a random walk in the energy space, e.g. in the second 
step of the cycle, system starts with its vacuum being, 
according to P^ . almost unpopulated and it has to re- 
lax to a state in which the probability of the ground state 
is almost 1. So, if every transition from one excited en- 
ergy level to another takes a finite time, one will need 
to perform in average n/2 transitions before visiting the 
ground state. Now note from (PT|) that the convergence 
of 77 to rjc is controlled by C'(l/ln[7i]): a smah step to- 
wards 77c will lead to a large increase in t^ nullifying the 
power C'(ln[77]/7T,) for n ^ 1; see (|20p . Hence for this type 
of relaxation the Carnot efficiency is not achievable at a 
finite power. 

The second relaxation scenario of the system with 
Hamiltonian (jll[ llSp is similar to the known combina- 
torial optimization problem: finding a marked item in an 
unstructured database [ll| of 77. -I- 1 items. This prob- 
lem is mapped to physics by associating each item to an 
eigenvector of a Hamiltonian. The marked item relates to 
the lowest energy level 0, while all other (excited) eigen- 
values of the Hamiltonian e are equal. The resulting sys- 
tem has unknown eigenvectors of the Hamiltonian, but 
known eigenenergies. Now the searching process can be 
organized as a relaxation of the system from an excited 
state to a low-temperature equilibrium state. This state 
is dominated by the ground level due to a large e. Once 
the relaxation is over the unknown itern (eigenvector) can 
be revealed by measuring the energy [15| . 

Within classical physics the minimal search time of this 
problem scales as 0{n) for n :^ 1 [ll|. It is thus not much 
better than going over all possible candidates for the so- 
lution, a typical situation of a computationally complex 
problem. For quantum algorithms (Grover's search) the 
search time scales as 0{y/n) [Uj. This is smaller than 
the classical estimate 0{n), but it still not suitable for 
our purposes, since it nullifies the power for In [77] ^ 1. 



Sub-optimal engine. Within the second (realistic) 
relaxation scenario, we shall modify the optimal engine so 
that the power is finite, but both the work and efiiciency 
are still large. We are guided by the analogy between 
the relaxation of the Hamiltonian ([TTl [T5|) under the sec- 
ond scenario and the Levinthal's paradox from protein 
physics [l^. In fact, pTl [Tsj) is the simplest model em- 
ployed for illustrating the paradox [13[ . Here the ground 
state refers to the unique folded (native) state. To ensure 
its stability, it must be separated by a large gap from ex- 
cited (free) energy levels. The essence of the paradox is 
that assuming the existence of many equivalent excited 
(unfolded) states, the relaxation time to the native state 
will be unrealistically long. Recall that the states pi and 
P2 of the optimal engine refer respectively to unfolded 
and folded states of the protein model; see after (HI]). 

The resolution of the paradox is to be sought via re- 
solving the degeneracy of excited levels: if there are en- 
ergy differences, some (unfavorable) transitions will not 



be made shortening the relaxation time J13l |. In resolv- 
ing the energy degeneracy we follow the basic model pro- 
posed for resolving the paradox [12|. 

The model has A^ ^ 1 degrees of freedom {(7^}"^]^; each 
one can be in C+ 1 states: 0-^ = 0, ..., (. Whenever Ci = 
for all i's, the model protein is in folded (ground) state 
with energy zero [ij, \lM- The ground state has zero 
energy. Excited states with s > 1 have energy e + Ss, 
where e > and s is the number of (misfolded) degrees 
of freedom with ai ^ 0. S > is the parameter that 
(partially) resolves the degeneracy of excited states; we 
revert to the previous, work-optimal model for (S — > 0. 
Thus for the Hamiltonians Hi and H2 we have 



i^a(.s) = [0, {fH+s^H} 



N 



1 2 



C°JV! 



(22) 



where each energy e^"! -f sj["l is degenerate gihy^Ui times. 
Given ([22]) , the functioning of the cycle consists of two 
isothermal and two thermally isolated pieces with sudden 
changes (jW,e'^l) ^ (-Jl^UI^l) ^ (j[il,e[il). Thus (HHIl 
remain true. Recall that they determine the generalized 
Carnot cycle that is optimal with respect to the thermally 
isolated pieces. Below we shall also assume 



/3i<5W=/32.5[2l, 



(23) 



because this makes the sub-optimal engine structurally 
very similar to the optimal one. Now the work W — 
Wi + W2 is calculated from (H H H [H [23l) : 

r.^ m(Ae+i|f#^)(.i~^2) ,^^. 

W[vi,V2]K] = ^-— i-- ^, 24 

(1 + mvi){l + m.V2) 

e[2l=Tiln[l/t;i]-r2ln[l/t;2], (25) 



Ae 



.[1] 



AS = (Jl^l - (S[21 = (Ti - T2) hi[C/K]. 



(26) 



where K = ^e 



',3i5W 



m 



(1 -I- if) , and where v^ 



-fS.e^' 



TABLE I: Parameters of the sub-optimal engine — work W, 
efficiency ry and the cycle time 2Tr — versus those of the opti- 
mal one; see (I24l428p . Wmax is the maximal work extracted for 



K = (e 



-ftilU 



at a vanishing power; see p6l I17p . A^ = 140, 



C = 4, Ti = 1, T2 = 1/2. Carnot efficiency r/c = 1/2. The 
Curzon-Ahlborn efficiency for this case is 0.585877c. 
Pi°' = [H-(H-is:)^e-''°''°']-^ (for Q = 1,2) are the ground- 
state probabilities of pa oc e~'^°^°; see (|22p . 



K 


Tr 


W/Wn... 


W 


v/vc 


P^^' 


[2] 


0.1 


4.45 X 10~^ s 


0.2267 


23.52 


0.8751 


0.0392 


0.9808 


0.2 


4.35 s 


0.3884 


40.3 


0.9110 


0.0237 


0.9883 


0.24 


357 s 


0.4393 


45.58 


0.9181 


0.0210 


0.9896 



Note the analogy between P^ and ([23]), with m being 
an analogue of n; they get equal for J — > 0. Note that in 
we neglected factor 0( — ) assuming that m^ 1. 
Likewise, we get for the efficiency: 



1 ^2 
77 = 1 

' Ti 



Ini 



NKln{v/K) 
l + K 



In^ 



NKln(v/K) 
l + K 



(27) 



For this model [12| assumes a local Markov relaxation 
dynamics, where each degree of freedom makes a tran- 
sition (Ji — > (Ti ± 1 in 10^^ seconds; this value is chosen 
conventionally to fit experimental magnitudes for the el- 
ementary dynamic step |12| . The model has a single re- 
laxation time (l2| that is easily reproduced in the gen- 
eral master-equation framework [section III of the Sup- 
plementary Material]: 



T, = 10~^(1 + KyyiNK) seconds. 



(28) 



(a = 1, 2) are determined from maximizing (|24l) . 



where the factor N is due to the A^-times degenerate first 
excited level. 

For (5["1 ^- (a = 1, 2), where the excited energy levels 
become degenerate, Tr ex (1 -|-C)^ scales linearly over the 
number of energy levels, as expected. When S^"' are not 
zero, Ti- can be of order of 1 second for N ^ 100, because 
l+K is close to 1. However, for the macroscopic situation 
[N ~ lO^'^) Tr is still huge. In this sense, the model is 
incoinplete, but still useful for analyzing the mesoscopic 
situation TV ~ 100 that is relevant for the protein folding 
problem [13[. 

Table I illustrates the characteristics of the sub-optimal 
engine and compares them with those of the optimal one. 
It is seen that a reasonable cycle times can coexist with 
a finite fraction (^ 40%) of the maximal work and with 
large efficiencies (^ 90%) of the maximal value that is 
definitely larger than e.g. the Curzon-Ahlborn value). 
Hence, albeit within the second (realistic) scenario it is 
impossible to approach the maximal efficiency as close 
as desired, reasonably large efficiencies at a finite (or 



even large) power are still possible. These results re- 
semble the power-efficiency trade-off (see [10|), but they 
are more complicated, since they involve work, efficiency 
and power. 

Summary. The generalized Carnot cycle allows to 
understand why the Carnot efficiency at a finite power is 
not realized in practice: it relates to the computational 
complexity that enforces very large cycle times. As we 
argued, this does not mean that finite-power and Carnot 
efficient engines are impossible to built in principle; it is 
that they are not met under realistic conditions. Even 
then, it is useful to study the optimal-work-efficiency 
engines: though they have an unrealistically long cy- 
cle time, their straightforward improvements do produce 
models working at reasonably large efficiency (and finite 
power). These sub-optimal models coincide with those 
of good folded proteins. We shall return elsewhere to 
biophysical implications of this fact. 

K.V.H. is supported by the Spanish project FIS2010- 
14830. 
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Supplementary material 

This material consists of four sections. Section I relates 
together two definitions of work. Section II estimates the 
relaxation time of the optimal engine within the master- 
equation framework (first scenario of relaxation) . Section 

III does the same for the sub-optimal engine. Section 

IV discusses an example of the first relaxation scenario 
and shows that the relaxation time of the optimal en- 
gine within this scenario is much shorter than the best 
estimate within the second scenario. 



I. CLARIFICATION OF THE CONCEPT OF 
WORK AS USED IN OUR SITUATION. 

Let a system interact with a source of work only (ther- 
mally isolated process). This means that the system 
Hamiltonian H[a{t)] is a function of a (classical) param- 
eter a{t). The work done on the system per unit of time 
equals to "force" times "displacement" and averaged over 
state of the system, as represented by a time-dependent 
density matrix pit): 



AW daU) fdH,^ 



(29) 



The same formula applies in the classical situation, where 
tr means integration over the phase-space (the space of 
coordinates and momenta), while p becomes the phase- 
space probability density. Now one can use the equations 
of motion for the density matrix, ih^ = H[oi{t)\p{t) — 
p{t)H[a{t)\ (in the classical situation this becomes the 
Liouville equation), to show from (|29p that the total work 
equals to the change of average energy 



W(t) 



At 



AW 
At 



tr(iJ[Q(r)]p(r)-i/[a(0)]p(0)). (30) 



II. RELAXATION TIME OF THE OPTIMAL 
ENGINE VIA MASTER EQUATION 

Consider a system with n ^ 1 degenerate (excited) levels 
with energy e > 0, and a single ground state with energy 
0. The relaxation time of this system can be estimated 
via the master equation. Let {pi}"^o be the probability 
of energy levels. The master equation reads 



Po 



En 
UlOiPi 
i— 1 



PO 



e: 



WiO, 



(31) 



where woi is rate of the transition i — >■ 0. Since all energy 
levels besides the lowest one have the same energy e, both 
woi = W(ii and Wio = wio do not depend on i. The rates 
relate to each other via the detailed balance conditions 
(which reflects the fact that the bath is in equilibrium at 
temperature 1//3): 



woie 



Wio. 



(32) 



In addition, since all the excited energy levels are equiv- 
alent, woi and uiio have an overall dependence 0{ — ): 



woi = win/n, wio = wio/n, 



(33) 



where i&io and woi do not depend on n. The origin of p3|) 
is that (apart of the energy difference between tlie ground 
state and an excited state) the transition takes place with 
a equal probabihty to (from) any excited energy level. 
We end up from (|31J433p with 



po = —W{nPo[ne '^^ + 1] + u;oi 



(34) 



meaning that the relaxation time Tr for poit) to converge 
to its equilibrium value 



,[<=q] 



pr = l/(l + ne-^') 



reads 



1/rr = ti)oi [ne ' +\]/n. 



(35) 



(36) 



We are interested in the situation, where the equilib- 
rium probability (|35p of the ground-state is sizable. More 
precisely in our situation e~^'^ = C'(ln[n]/n) or e~^^ = 



0{ 



). At any rate, the relaxation time n- in 



n ln[n] 

roughly 0{n) 



III. RELAXATION TIME OF SUB-OPTIMAL 

ENGINES VIA ADIABATIC APPROXIMATION 

OF THE MASTER EQUATION 

We return to (|3H) . but we do not assume anymore that 
the excited levels ei, ..., £„ have the same energy. But we 
still assume that the gap £i > between the zero-energy 
ground state and the first excited state is the largest en- 
ergy parameter in the system, because we want to ensure 
that the equilibrium ground-state probability is close to 
1. Thus we may apply to (|3ip the adiabatic approxima- 
tion assuming that on those times where po changes, the 
excited-state probabilities pi already equilibrated: 



p,(t) = (l-po(i))e- 



'e: 



-/3e, 



(37) 



Hence the relaxation time Tr of po deduced from pil I37p 
reads 






Woi 



(38) 



where the employed the detailed balance condition. 
Eq. psp reproduces the relaxation time (|28p of the 
Zwanzig model [Q, if we employ the energy spectrum 
(|22p and note that the transition probabilities from the 
first excited (("A'^-degenerate) energy level are constants, 
wok = 10® for k = 1,...,(^N, while no transitions (to 
the ground state) is possible from other excited states: 
woi — for I = i^N + Ij ..., n. Here 10® is the charac- 
teristic microscopic scale [12] . Putting these into (|38p we 
get that the first term in the right-hand-side of (|38p re- 
produces (|28p . The second term is neglegible, if e~^^^ is 
sufficiently small. 

Note that according to psp . the relaxation time of the 
Zwanzig model is still unacceptably large, if A^ 3> 1. It is 
not difficult to get rid of this restriction, though it is not 
useful from the viewpoint of the heat engine functioning, 
where a large parameter is needed. 



IV. AN EXAMPLE OF THE FIRST 
RELAXATION SCENARIO 

System-bath interaction 

Consider a system E with n excited energy levels and one 
lowest energy (ground state) whose energy we set to zero. 
All n excited levels have the same energy e. 
The initial density matrix of E is Gibbsian at tempera- 
ture To = l//3o: 



p oc e ' " = rVu H Vc 

n 

1 



(39) 

■ 1 + ne-^o. (40) 

where Vo = |0)(0| is the projector on the ground state, 
and Pe is the projectors on the n-dimensional eigen-space 
of p with eigenvalue e. 

Now E interacts with an external thermal bath at tem- 
perature T — 1/P, so that the density matrix of E con- 
verges in time to peq oc e~^^ . We shall design a concrete 
model for this interaction and estimate the relaxation 
time. 

We assume that the bath consists of a large number 
of independent particles prepared in identical (thermal 
states). E interacts with one particle, then with the sec- 
ond one etc. Since the particles are independent, it will 
suffice to consider the interaction of E with the first par- 
ticle B only. 

We assume that the bath particle B has (among other 
energies) energy levels E and E + e. The degeneracies 
of these levels are ue and ns+e, respectively. The initial 
(before interacting with E) equilibrium density matrix of 
B reads 



-I3Hb 



Ve 



/Z = (J + rE'RE-\- r^+E H 



-I3E 



/z, 



c+E, 

E-HE 
riEe "^ , 
E 



(41) 
(42) 



where te and r_B+£ are the Boltzmann weights for the 
energy levels E and E + e, respectively, the summation 
in (|42p is taken over all energy levels of B. H_b and Tle+E 
are the projectors on the corresponding sub-spaces, 

trllE = nE, tiUs+E = nE+e, (43) 



and where a in (|4ip is the remainder of a. 
It is assumed that the unitary operator V responsible 
for the interaction operates within the sub-space with 
the projector Ve €) Hb -(- Po ® H_b+e (this sub-space has 
energy E + e), i.e.. 



[v,re(»iiE + ro(»'nE+e] =o. 



(44) 



On the remainder of the overall Hilbert space (of E-l-B) 
V acts as unit operator. Thus, V commutes with the 
Hamiltonian of E-l-B. Hence no additional energy (work) 
is nedeed for switching the E-B interaction on and off. In 
that respect V resembles the weak-coupling, though by 
itself it need not be weak, i.e. it need not be smaller than 
the Hamiltonian of E-f B. 
Then the post-interaction density matrix p' of E reads 

p' = tTsVpCSaV^ 

( 1- 

= p — r TE+e — TE 



(45) 



[riE+eVo - trfiV-po ® IlE+eV^]. 



It will suffice to keep track of the lowest energy-level oc- 
cupation (0|p'|0) = r' only: 



r' -r = ~A[r~T,^], r,^^ -—-—^^, (46) 

A=^^ \nE+, - (01 ( treVPo ® n^+^Vt )10)] , (47) 

r-cqn L J 

where r^q is the equilibrium value of r. Using (|42|l one 
can show that A < ^max < 1: after the interaction E 
gets closer to its equilibrium state; see (|46p . 
Now (|45|l serves as the initial state of E for a similar 
interaction with the second bath particle that initially 
has the same state a as in (|41|l . We get for all subsequent 
interactions fwe revert from (1481) to (1461) for m — 1]: 



M] 



= {l- A)"" [r - 



(48) 



It is seen that (|48|l predicts exponential (with respect to 
the number of collisions) relaxation towards the equilib- 
rium value roq of r. The approach to equilibrium is gov- 
erned by the factor (1 — ^4)™ meaning that when |A| ^ 1 
the effective number of interactions after which the equi- 
librium is established (which is proportional to the relax- 
ation time) equals to — l/[ln(l — A)]. 



Minimization of the relaxation time 

Since we are interested in possibly shorter relaxation 
time, we need to maximize A over the unitary V [under 
condition (|44|l ]. To do that, we first write Po<8>n_E+£ in a 
conveniently chosen matrix form in the energy eigenbasis 

PolSlTlE+e = (49) 

diag(...,l,...,l,...i...,0,...,0,...j---|...,0,...,0,...),(50) 



^R^s elements 



n sections 



where sections correspond to eigenvectors of H , and el- 
ements in sections run over the ones of Hb- Values are 
shown only for the subspace given by Po^Ti-E+e+Pe'SiTlE- 
The first section (denoted as |...|) there are ue+e unities. 
Then come n identical sections, each one contains he 
zeroes. 

To maximize A we need to minimize (0| {trsVPo <8 
TIe+eV'' )\0) over all possible Vs living in Pq <8) ^e+e + 
Pe^TIe- One can show that the optimal unitary amounts 
to a permutation of the eigenvalues (|50p (this fact can 
be shown similarly to the derivation presented in (f9l ll3p 
of the main text). Now we note that the trace over tr^ 
amounts to summing up elements in each section. So the 
element (0| {trsVPo^'n-E+EV^ )|0) wiU be the sum of the 
elements of the first section in the permuted diagonal in 
(|49|) . The optimal permutation will thus be the one which 
takes out of the first section as much unities as possi- 
ble. Therefore, if ue+e (the number of unities) is < than 
UEn — the number of zeroes — then it is possible to move 
all unities, making miny {{0| (trsVPo ® TIe+eV^ )|0)} ~ 
0. Otherwise, the latter quantity will be ue+e — riETi, 
leading us to the following formula: 



^max — 



te min [ ue+e , ueti ] 



(51) 



Relaxation time for realistic bath 

Let us work out (|51|) for a realistic example of the bath. 
We assume that the bath particle amounts to L 3> 1 
independent two-level systems. Each such system has 
energies and C,. Thus the bath particle has energies 
0, C,, 2C,, ..., LC,. Each bath energy level E is degenerate 



riE = 



L\ 



{E/0\{L-EK)\ 
times. Provided that 



Ue+e > UEn, 



we obtain 



Sex nEC 



-I3E 



Amax = (1 + ne' 
Since we want a larger Amax, we take 

E = L/{e'"^ + l). 



(52) 



(53) 



(54) 



(55) 



Hence using the Stirling's formula L\ ~ \J2-kLL^ e ^ and 
noting that Z = {1 + e"''^)^ we get from ((SJ) 



(1 -f ne 



-;9.^ 



H-e-« 



(56) 



We work out (|53p via the Stirling's formula and obtain 
from dia [53 



ln[n] 
L 



h2 



1 



g/3C + l 



</l2[ 



e0i + l L(' 



where h2[x] = — a;lna:: — (1 — a;) ln(l — x). 
Let us specify e as [^ is a parameter] 

e = ^Tln[n], ^ > 1. 

Putting this into (|57|l we get 



ln[n] 



+ /i2[- 



Pi + l 



]<h2[ 



+ ii .^'^t" 



e^C + 1 PC L 



(57) 



(58) 



(59) 



Provided that /i > 1, (|59p can be satisfied for sufficiently 
small (but finite) -2^^ and sufficiently large (" > 0. Note 
that ()59p never holds for fJ, < 1. 

Returning to (|56p we see from (|58p that for ln[n] = 
0{L) > 1 



Oi- 



1 



(60) 



which means that the relaxation time scales as 

o(\/H^). 

Consider now the opposite [to (|53p ] case 



where 

^niax 

Choosing 



Ue+e < nEU, 



(1 + r::^)^^^ 



-/3(B+£) 



-Pe' 



iJ + e = L/(e'^^ + l), 



(61) 



(62) 



(63) 



we get [cf. (f56)) ] 



Relations with the main text 



1 N /1 + e-'"; 



If we specify [cf. (f55)) ] 

e = uT\n[n], v < I, 
then ([SI]) reads [cf. ((59l) ] 



, r 1 1 ln[nl , , 1 



u ln[ri 



'e/3« + l' i ^e-sc + l /3(; L 



(64) 



(65) 



(66) 



This relation holds for v < 1, sufficiently small (but fi- 
nite) isM and sufficiently large C > 0. 
Hence from (|66l I62|l we return to the same conclusion 

(Eni. 



In the main text — see in particular (|16l I19p — we stud- 
ied the relaxation of the optimal engine E that has the 
energy spectrum described at beginning of section IV. 
More specifically, in the main text we needed two differ- 
ent relaxation scenario: E with energy gap e'^' ~ Ti ln[n] 
relaxes on a thermal bath at temperature T2 (where 
Ti > T2), and conversely E with energy gap e'^' ~ T2 ln[n] 
relaxes on a thermal bath at temperature Ti ; see (|19p of 
the main text in this context. 

Now the first case corresponds to (|58l 1591 I60p . while 
the second case to (|65l 1661 I60p . In both cases we get 
that for a realistic thermal thermal bath (but with tuned 
system-bath interactions) the relajcation time amounts to 
0(\/InR); see ^. 



